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second  probe  mode  with  strong  atom-photon  interactions.  A  new  hybrid  trap 
is  analyzed  that  combines  optical  and  Casimir-Polder  forces  to  form  stable 
traps  for  neutral  atoms  in  dielectric  nanostructures.  By  suitable  design  of  the 
band  structure,  the  atomic  spontaneous  emission  rate  into  the  probe  mode  can 
exceed  the  rate  into  all  other  modes  by  more  than  tenfold.  The  unprecedented 
single-atom  reflectivity  ro  >  0.9  for  the  guided  probe  field  should  enable  diverse 
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1.  Introduction 

New  opportunities  in  atomic,  molecular  and  optical  (AMO)  physics  and  quantum  information 
science  emerge  from  the  capability  to  achieve  strong  radiative  interactions  between  single 
atoms  and  the  fields  of  nanoscopic  optical  waveguides  and  resonators  [1].  For  example, 
strong  atom-photon  interactions  in  lithographic  structures  [2-6]  could  be  used  to  create 
quantum  optical  circuits  with  long-range  atom-atom  interactions  mediated  by  single  photons 
[7,  8].  Moreover,  linear  arrays  of  atoms  radiatively  coupled  to  nanophotonic  waveguides 
exhibit  a  wide  range  of  remarkable  phenomena,  including  coherent  transport  of  atomic 
emission  [9-11],  guided  superradiance  and  polaritons  [12-14],  as  well  as  highly  reflecting 
atomic  mirrors  [15,  16].  The  interplay  of  atomic  emission  into  the  waveguide  and  photon- 
mediated  forces  can  lead  to  self-organization  of  atoms  into  exotic  spatial  configurations  along 
the  waveguide  [17,  18]. 

A  long-standing  obstacle  to  this  scientific  frontier  is  the  challenge  of  trapping  atoms  in 
vacuum  near  dielectric  surfaces  (~100  nm)  while  at  the  same  time  achieving  strong  interactions 
between  one  atom  and  photon.  A  far-off  resonance  dipole-force  trap  (FORT)  [19]  can  provide 
atomic  localization  by  using  modes  of  the  dielectric  for  optical  trapping  [20-22]  and  has  been 
used  to  trap  cold  atoms  within  hollow-core  optical  fibers  [23-26]  and  external  to  fiber-taper 
waveguides  [27-29]. 

Motivated  by  these  advances,  in  this  paper  we  present  principles  for  the  design  of 
optical  traps  and  strong  atom-photon  interactions  in  one-dimensional  (ID)  photonic  crystal 
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Figure  1.  (a)  Schematic  for  the  single  nanobeam  structure  with  dimensions 
(a,  w,  t,  hx,  hy )  =  (367,  845,  825,  246,  745)  nm.  (b)  Band  diagram  for  the  single 
nanobeam  in  (a)  showing  only  bands  with  even  vector  symmetry  about  the  y 
and  z  symmetry  planes.  The  trapping  and  probing  bands  are  shown  as  thicker 
lines,  with  the  trap  ojj/2tz  (probe  oja/2tt)  frequency  as  a  blue  (red)  dashed  line, 
(c)  Field  intensity  of  the  blue  trapping  mode  and  (d)  field  amplitude  of  the  probe 
mode  in  the  center  plane  z  =  0  for  the  single  nanobeam  in  (a).  Green  spheres 
mark  the  locations  of  minima  of  the  trapping  potential,  (e)  Schematic  for  the 
double  nanobeam  structure  with  (a,  w,  t,  d,  g )  =  (335,  335,  200,  116,  250)  nm. 
(f)  Band  diagram  for  the  double  nanobeams  in  (e)  displaying  only  modes  of  even 
vector  symmetry  in  z.  The  proximity  of  the  two  nanobeams  results  in  a  band 
structure  composed  of  even  (green)  and  odd  (magenta)  superpositions  of  single 
nanobeam  modes.  We  focus  on  the  even  parity  supermodes  due  to  their  large  field 
amplitude  in  the  gap.  (g)  Field  intensity  of  the  blue  trapping  mode  and  (h)  field 
amplitude  of  the  probe  mode  in  the  center  plane  z  =  0  of  the  double  nanobeams 
in  (e).  The  black  diamonds  in  (b)  (f)  mark  resonances  for  finite  structures  of 
81  unit  cells  from  figures  4  and  5). 


waveguides.  Our  choice  of  ID  waveguides  is  based  on  their  simple  photonic  band  structures 
and  their  compatibility  with  current  laser  cooling  and  trapping  technology  [27-29]. 

We  analyze  the  trap  potential  Ut0 t(r)  due  to  light-shifts  from  a  FORT  [30-32]  together 
with  Casimir-Polder  (CP)  interactions  with  the  dielectric  [33-36]  (figures  1(a)  and  (e)). 
Despite  the  proximity  of  the  surfaces,  stable  potentials  Ut0 t(r)  are  achieved  for  modest  optical 
intensities  (~5  mW  /xm-2)  for  blue-detuned  FORTs  operated  at  a  ‘magic’  wavelength  for 
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the  D2  line  of  atomic  cesium  [19].  A  new  possibility  for  trapping  is  also  identified  for 
which  vacuum  forces  from  CP  interactions  are  exploited  to  close  the  trap  perpendicular  to 
the  plane  of  structure,  which  would  otherwise  be  unstable  with  either  the  FORT  or  the  CP 
potential  alone  [37], 

In  addition  to  the  waveguide  trapping  properties,  strong  near-resonance  atom-photon 
interactions  of  trapped  atoms  are  found  to  arise  for  waveguides  with  properly  tailored  band 
structure  [38-42],  For  practically  realizable  structures,  we  find  yw/y'  >  10,  where  )/|D  is  the 
atomic  decay  rate  into  the  (guided)  probe  mode  and  y'  the  rate  into  all  other  modes.  One 
atom  trapped  within  the  structure  could  thereby  attenuate  a  resonant  probe  with  transmission 
|l-r0|2<  10“2[9,  16]. 


2.  Design  principles  for  the  one-dimensional  photonic  crystal  waveguides 

As  illustrated  in  figure  1,  we  focus  on  two  of  the  simplest  quasi- ID  photonic  crystal  geometries. 
The  first  waveguide  consists  of  a  single  silicon-nitride  nanobeam  (refractive  index  n  =  2)  with 
a  ID  array  of  filleted  rectangular  holes  along  the  propagation  direction;  atoms  are  trapped  in  the 
centers  of  the  holes  (figure  1(a)).  The  second  waveguide  consists  of  two  parallel  silicon  nitride 
nanobeams,  each  with  a  periodic  array  of  circular  holes,  with  atoms  trapped  in  the  gap  between 
the  beams  (figure  1(e)). 

The  design  of  a  lD-photonic  crystal  waveguide  with  distinct  modes  for  optical  trapping 
and  strong  atom-photon  interactions  is  constrained  by  the  region  of  the  optical  band  structure 
containing  a  continuum  of  unguided  optical  modes  (i.e.  the  light  cone  indicated  in  gray 
in  figures  1(b)  and  (f)).  Modes  within  the  light  cone  can  still  have  large  amplitude  in  the 
structure  but  radiate  energy  into  the  surrounding  vacuum  leading  to  unacceptable  loss.  The 
top  of  the  vacuum  light  line  is  at  the  Brillouin  zone  boundary  (A -point,  where  kxa  =  Jt), 
so  the  lattice  constant  a  is  constrained  by  a  <  A/2,  where  A.  is  the  smaller  of  the  (vacuum) 
wavelengths  for  trapping  and  probe  fields.  Here,  kx  is  the  Bloch  wavevector  along  the  waveguide 
axis  x. 

Once  a  is  fixed,  additional  guided  modes  can  be  ‘pulled’  below  the  light  line  by  increasing 
the  width  and  thickness  of  the  structure.  With  appropriate  modes  below  the  light  line  for  probing 
and  trapping,  the  spacing  of  these  modes  at  the  A -point  can  be  tuned  by  altering  the  size  of 
the  holes,  which  enables  the  probe  mode  to  be  resonant  with  the  frequency  coA  of  the  atomic 
transition  while  simultaneously  matching  the  optical  frequency  ojj  of  the  trap  mode  to  a  ‘magic’ 
frequency  for  the  atom  [19]. 

3.  Optical  trapping  and  surface  Casimir-Polder  potential  inside  the  ID-waveguides 

Within  this  general  context,  here  we  consider  only  blue-detuned  FORTs  for  which  the  trapping 
mode  has  an  intensity  minimum  at  the  trapping  site7.  Our  analysis  is  for  the  D2  line  of  atomic  Cs 
with  probe  wavelength  near  the  atomic  resonance  aa  =  852  nm  and  with  a  blue-detuned  FORT 
at  the  magic  wavelength  A.T  =  793  nm  [19].  Note  that  our  results  are  readily  transcribed  to  other 
atomic  transitions  by  way  of  the  scale  invariance  of  Maxwell’s  equations  [43]. 

The  photonic  crystals  are  assumed  to  be  suspended  in  vacuum  and  composed  of  SiN. 
Band  structures  are  calculated  using  the  MIT  photonic-bands  (MPB)  software  package  [44], 

7  We  have  also  designed  structures  for  trap  and  probe  modes  for  red-detuned  FORTs  with  an  atom  trapped  at  the 
intensity  maximum  in  a  unit  cell. 
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Figure  2.  Trapping  potentials  for  the  single  nanobeam  structure  in  figure  1(a) 
for  Cs  6S1/2,  F  =  4  level  and  at  =  793  nm.  (a)  CP  potential  Uc p(r)  and  (b)  total 
potential  f/tot(r)  =  Uc p(r)  +  t/F oRi(r)  in  the  central  z  =  0  plane.  Panels  (c)-(e) 
show  line  cuts  of  Uc p  (red  solid),  UF ort  (blue  dashed)  and  Utot  (blue  solid) 
along  the  (c)  x-,  (d)  y-  and  (e)  z-axis.  Average  trap  intensity  for  a  unit  cell  is 
4.9  mW  \i  m-2. 


Field  profiles  for  guided  modes  are  calculated  using  the  finite-element-method  simulations8. 
Results  for  the  single  and  double  nanobeam  structures  are  presented  in  figure  1 . 

With  suitable  guided  modes  for  trapping  in  hand,  we  have  developed  numerical  tools  for 
evaluating  the  FORT  and  CP  potentials  inside  the  waveguide,  and  hence  the  total  potential 
C/tot  =  U fort  +  Uc p-  The  adiabatic  potential  U fort  ( r)  is  readily  calculated  using  the  electric  field 
distribution  of  the  trap  mode,  E(r)  =  U/,r  (r)  e'kxX  [30-32];  U/,A  (r)  is  the  periodic  Bloch  wave 
function  at  kx . 

The  surface  potential  Uc p(ra)  is  determined  from  the  formalism  in  [35]  for  the  imaginary 
component  of  the  scattering  Green’s  tensor  Gsc(ra,  ra,  co)  from  Maxwell’s  equations  for  a  point 
dipole  at  the  atomic  location  ra  with  the  vacuum  contribution  (i.e.  no  dielectric  structure) 
subtracted.  We  evaluate  Gsc(ra,  ra,  co)  numerically  by  adapting  the  procedures  from  [36]  (see 
appendix  A). 

Figures  2  and  3  display  numerical  results  for  Ccp(r),  C/fortW  and  f/tot(r)  for  the  single 
and  double  nanobeams  for  kx  below  the  X -point.  The  calculations  are  for  the  6S1/2,  F  =  4 
hyperfine  ground  state  of  Cs  for  the  FORT  modes  indicated  in  figures  1(c)  and  (g)9.  For  these 

8  COMSOL  Multiphysics  (http://comsol.com/). 

9  Utot(r)  is  largely  insensitive  to  the  particular  atomic  Zeeman  sublevel  owing  to  negligible  vector  shifts  due  to  the 
high  degree  of  linear  polarization  of  the  trapping  modes  and  to  the  absence  of  ground  state  tensor  polarizability. 
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Figure  3.  Trapping  potentials  for  the  double  nanobeam  structure  in  figure  1(e) 
for  Cs  6Si/2,  F  =  4  level  and  at  =  793  nm.  (a)  CP  potential  Uc p(r)  and  (b)  total 
potential  Utot(r)  —  £/cp(r)  +  C/fortW  in  the  transverse  a  =  0  plane.  Panels 
(c)-(e)  show  line  cuts  of  Ucp  (red  solid),  C/fort  (blue  dashed)  and  Um  (blue 
solid)  along  the  (c)  x-,  (d)  y-  and  (e)  z-axis.  Average  trap  intensity  for  a  unit  cell 
is  3.5  mW  /x m-2. 


initial  calculations,  we  make  the  reasonable  assumption  for  SiN  that  the  dielectric  constant  € 
is  frequency  independent,  e(r,  co)  — >  e(r);  see  appendix  A.l  for  the  justification  by  numerical 
evaluation  of  £/Cp(r)  using  the  full  frequency  dependence  of  e(co)  for  SiN  at  several  locations. 

3.1.  Stable  optical  trapping  in  a  single  nanobeam 

Ut ot(r)  for  the  single  nanobeam  in  figure  2  reveals  that  modest  optical  intensity  is  sufficient 
to  overcome  the  attractive  CP  interactions  and  create  a  stable  potential  in  the  center  of  the 
vacuum  space  at  rmin  =  0  within  a  unit  cell.  An  atom  would  be  localized  at  distances  (dx ,  dy)  = 
(123,  373)  nm  from  the  walls  of  the  dielectric.  The  trap  oscillation  frequencies  for  a  Cs  atom 
would  be  (fx,  fy,  fz)  ~  (612,  180,  484)  kHz. 

3.2.  A  new  hybrid  trap  in  double  nanobeams:  utilizing  vacuum  Casimir-Polder  force  to  realize 
a  stable  trap  condition 

For  the  double  nanobeams  [45],  the  FORT  alone  is  insufficient  to  trap  the  atom,  as  the  mode  has 
a  (weak)  local  intensity  maximum  along  the  z  direction  that  repels  an  atom.  However,  the  CP 
potential  Uc p(r)  along  z  provides  the  force  necessary  to  overcome  the  repulsive  optical  force 
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and  to  form  a  stable  trap.  The  result  is  a  hybrid  optical-vacuum  trap  that  circumvents  the  ‘no-go’ 
theorem  for  vacuum  trapping  alone  [37], 

Potentials  for  our  hybrid  trap  are  illustrated  in  figure  3.  At  the  trap  minimum  rmm  =  0,  an 
atom  would  be  localized  at  distance  dy  —  125  nm  from  adjacent  surfaces  of  the  dielectric  beams. 
Oscillation  frequencies  for  a  Cs  atom  would  be  (fx,  fy,  fz )  ~  (1013,  390,  57)  kHz. 

4.  Achieving  strong  atom-photon  interactions  via  guided  modes  near  a  photonic 
band  edge 

As  concerns  strong  radiative  interactions,  our  structures  trap  an  atom  in  a  region  of  large 
amplitude  for  the  probe  field,  leading  to  small  mode  volume  per  unit  cell  (appendix  E). 
Atom-photon  interactions  can  be  further  enhanced  near  a  band  edge  [38-42],  where  the 
density  of  states  diverges  due  to  a  van  Hove  singularity.  To  quantify  the  radiative  coupling, 
we  determine  the  decay  rate  ylot  for  a  point  dipole  located  at  ra  =  0  for  a  structure 
with  N  unit  cells10.  Finite-difference-time-domain  (FDTD)  calculations  are  performed  to 
evaluate  the  classical  Green’s  tensor  G(ra,  ra,  oj)  and  thence  ytot  following  [34,  35,  46] 
(appendix  B). 

4.1.  Green’s  tensor  and  atomic  spontaneous  emission  inside  the  ID-waveguide 

Figures  4(a)  and  5(a)  display  the  diagonal  components  Im[G„(vd)]  of  the  Green’s  tensor  as 
functions  of  dipole  frequency  vd  —  coA/2n  and  relate  to  the  emission  rate  of  resonant  point 
dipoles  polarized  along  the  i  —x-,  y-  or  z-axis  for  the  single  and  double  nanobeams.  Firstly,  in 
figure  4(a),  Im[G„]  is  enhanced  along  the  x  (periodic)  direction  across  a  broad  frequency  range, 
and  is  suppressed  in  the  y  and  z  directions,  as  can  be  explained  by  the  orientation  of  the  induced 
array  of  image  dipoles  along  the  single  nanobeam. 

Secondly,  figures  4(a)  and  5(a)  display  a  series  of  resonant  peaks  due  to  strong  emission 
into  various  guided  modes.  In  the  region  near  the  Cs  D2  line  (i.e.  vd  ~  vA  =  oja/2tt  =  352  THz), 
we  find  peaks  in  Im[Gxx]  for  the  single  nanobeam  and  in  Im[G  vv]  for  the  double  nanobeams. 
These  peaks  are  due  to  emission  into  our  designated  probe  modes  for  the  respective  structures, 
where  for  the  single  (double)  nanobeam(s),  the  probe  mode  is  principally  polarized  along  the 
x-(y-)  axis.  Each  peak  is  from  a  discrete  set  of  propagation  constants  ki")  ~  nn/aN  imposed 
by  the  boundary  conditions  for  the  finite  structures.  Here,  N  is  the  total  number  of  cells  in  the 
single  (double)  beam,  and  n  ^  N  is  an  even  (odd)  integer.  We  find  excellent  agreement  between 
the  frequencies  of  these  resonances  and  the  band  diagram  of  the  probe  mode  (‘diamonds’  in 
figures  1(b)  and  (f))  for  various  (n,  N). 

The  peaks  become  larger  and  narrower  as  k(xn)  approaches  the  A -point,  owing  to  the 
diminishing  group  velocity  [40,  42].  Beyond  the  A -point,  the  probe  resonances  disappear, 
leaving  a  broad  background  corresponding  to  coupling  into  lossy  (radiation)  modes. 

On  an  expanded  frequency  scale  around  vd  ~  va,  figures  4(b)  and  5(b)  show  calculated 
atomic  decay  rates  ytot  for  the  6P3/2,  F'  =  5  — >  6Si/2,  F  =  4  transition  in  atomic  Cs  (appendix  B). 
When  the  atomic  dipole  is  aligned  along  the  principal  polarization  of  the  designated  probe 
mode  (x — single  beam;  y — double  beam),  the  emission  rate  y1D  into  the  probe  mode  is  strongly 

10  A  finite  number  of  unit  cells  is  chosen  to  better  represent  an  actual  structure.  Since  our  focus  is  not  on  specific 
boundaries,  the  structure  is  terminated  by  vacuum. 
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vd  (THz) 


vd  (THz) 


Figure  4.  Green’s  tensor  and  total  atomic  decay  rate  ytot  versus  source  dipole 
frequency  vd  for  the  single  nanobeam  at  r  =  0.  Panel  (a)  shows  the  diagonal 
components  of  the  Green’s  tensor,  Im [Gxx]  (solid  black),  Im[GvvJ  (red)  and 
Im[Gzz]  (blue),  normalized  to  the  free  space  value  Im[G0]  (dashed  line).  The 
number  of  unit  cells  is  N  —  81.  Cesium  D2-line  frequency  vA  =  352  THz  is 
centered  in  the  shaded  area.  The  vertical  dotted  line  marks  the  light  line,  beyond 
which  all  decay  channels  are  lossy.  Panels  (b)  and  (c)  show  ytot  (black  curves), 
normalized  to  the  free  space  value  yo  (dotted  line),  in  the  frequency  range  marked 
by  the  shaded  area  in  (a).  The  solid  (dashed)  curve  is  evaluated  using  81(61)  unit 
cells.  The  atomic  spin  is  aligned  to  the  v-axis,  with  the  spin  projection  quantum 
number  (b)  mF  =  0  and  (c)  mF  =  5. 


enhanced  at  frequencies  corresponding  to  k[n>  near  the  X -point.  Specifically,  for  vd  =  vA  large 
enhancements  in  yjD  occur  for  the  initial  excited  state  6P3/2,  F'  =  5,  mF  —  0,  while  yiD  is 
suppressed  for  the  initial  state  m'F  =  5.  This  is  because  the  probe  mode  predominantly  supports 
n -polarization  and  hence  A mF  =  0.  Coupling  between  states  with  AmF  #  0  is  small.  Of  course, 
additional  guided  modes  can  contribute  to  yto t,  as  is  evidenced  for  the  mF  =  5  state  due  to  field 
polarizations  perpendicular  to  the  atomic  spin,  such  as  z  (Jc)  for  the  single  (double)  beam(s)  in 
figure  4(a)  (figure  5(a)). 

4.2.  Enhanced  probe-mode  coupling  rate  y\D  and  the  resulting  strong  atom-photon  interaction 

From  y’tot(fd)  and  an  analytic  model  of  coupling  to  the  guided-mode  near  the  X -point,  we 
estimate  the  contributions  of  ym  and  y'  to  ytot  =  yiD  +  y'  near  the  largest  resonances  in 
figures  4(b)  and  5(b)  (appendix  D.l).  For  mF  =  0  and  N  =  81,  we  find  that  ym/Yo  —  15 
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Figure  5.  As  in  figure  4  now  for  the  double  nanobeams  at  r  =  0.  Panel  (a)  shows 
the  diagonal  components  of  the  Green’s  tensor,  Im[Gxx]  (solid  black),  Im[Gvv] 
(red),  and  Im[Gzz]  (blue),  normalized  to  Im[Go]  for  free  space  (dashed  line). 
N  =  81  unit  cells.  Cesium  D2-line  frequency  va  =  352 THz  is  centered  in  the 
shaded  area.  The  vertical  dotted  line  marks  the  light  line,  beyond  which  all  decay 
channels  are  lossy.  Panels  (b)  and  (c)  show  ytot  (red  curves),  normalized  to  the 
free  space  value  y0  (dotted  line),  in  the  frequency  range  marked  by  the  shaded 
area  in  (a).  The  solid  (dashed)  curve  is  evaluated  using  81(61)  unit  cells.  The 
atomic  spin  is  aligned  to  the  y-axis,  with  the  spin  projection  quantum  number 
(b)  and  (c)  m'F  —  5. 


and  y'/Xo  —  1.2  for  the  single  nanobeam,  while  y1D/y0  ~  21  and  y'/y0  ~  1.0  for  the  double 
nanobeams11.  Here,  Yq/2tt  —  5.2  MHz,  the  free-space  Cs  decay  rate. 

The  ratios  yw /  ytot  and  yiD  /  y '  serve  as  metrics  for  the  strength  of  atom-photon  interactions 
for  our  ID  photonic  crystals.  For  example,  the  resonant  reflectivity  ro  of  a  trapped  atom  for 
the  probe  field  should  scale  as  r0  =  Xm/Xtot  [9,  16],  which  for  the  double  nanobeams  leads 
to  r0  ~  0.95.  For  a  cavity  quantum  electrodynamics  (QED)  system  with  one  ‘impurity’  atom 
surrounded  by  NA  ‘mirror’  atoms  along  a  lD-lattice  [16],  the  ratio  of  the  coherent  coupling 
rate  gi  =  v/NaXid/2  to  the  effective  dissipative  rate  y'  would  exceed  unity  even  for  NA  =  1 
atom.  For  conventional  cavity  QED,  we  estimate  a  1 -photon  Rabi  frequency  ~27T  x  2  GHz  for 
figure  5(b)  with  N  —  81. 

11  For  comparison,  ym/yo  —  0.06  is  the  state  of  the  art  for  a  Cs  atom  trapped  near  an  SiC>2  nanofiber  [29]. 
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5.  Conclusion 

Certainly  there  are  challenges  in  the  implementation  of  our  designs,  including  atom  loading 
into  the  small  trap  volumes  and  light  scattering  from  device  imperfections.  The  latter  can 
reduce  the  expected  coupling  rate  ym  and  introduce  trap  imperfections.  The  former  challenge 
may  be  resolved  either  by  loading  atoms  continuously  from  an  atomic  reservoir  surrounding 
the  nanostructure  [27-29]  to  improve  the  loading  efficiency,  or  by  using  a  moving  optical 
lattice  to  transfer  atoms  precisely  into  the  trapping  sites.  We  are  working  to  address  these 
issues  by  numerical  simulation,  device  fabrication  and  cold-atom  experiments  with  nanoscopic 
structures.  Our  efforts  are  motivated  by  the  prediction  yw/Ym  J  0.9  in  figures  4  and  5,  which 
is  unprecedented  in  AMO  physics  and  which  could  create  new  scientific  opportunities  (e.g. 
quantum  many-body  physics  for  ID  atom  chains  with  photon-mediated  interactions,  and  high- 
precision  studies  of  vacuum  forces).  Moreover,  our  double  nanobeam  structure  provides  proof- 
of-principle  for  a  promising  new  concept  that  combines  optical  and  vacuum  forces  to  form  stable 
traps  for  neutral  atoms  in  dielectric  nanostructures. 

Our  designs  of  ID  periodic  structures  represent  a  basis  for  extensions  to  more  complicated 
geometries,  including  two-dimensional  (2D)  photonic  crystals  with  trapped  atoms.  A  step  in  this 
direction  is  with  structures  that  are  still  ID  in  nature  (e.g.  a  line  of  defect  patterns  embedded  in 
a  2D  photonic  band  gap  crystal).  Such  structures  should  offer  larger  values  for  ym  and  further 
suppressions  of  emission  rate  y'  to  other  modes  (including  the  vacuum). 
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Appendix  A.  Calculation  of  Casimir-Polder  potentials 


The  CP  potential  t/cp(ra)  is  calculated  from  the  following  integral  [35]: 

r  pOO 

I  d  coco2 

Jo 


tlULn 

CCp(ra)  =  —  Im 


2n 


(A.l) 


where  Tr[.]  denotes  the  trace,  ar°  is  the  dynamic  polarizability  tensor  of  ground-state  cesium 
atom,  and  Gsc(ra,  ra,  co)  =  G(ra,  ra,  co)  —  Go(ra,  ra,  co)  is  the  scattering  Green’s  tensor,  that  is, 
the  Green’s  tensor  G  subtracted  by  the  vacuum  contribution  G0  evaluated  at  atomic  location  ra; 
2nh  is  Planck’s  constant,  and  //0  is  the  vacuum  permeability.  The  Green’s  tensor  is  the  solution 
to  the  Maxwell  equation  [V  x  V  x  —  ^-e(r,  oj)\G(r,  r',  co)  =  I6(3)(r  —  r'),  corresponding  to  the 
electric  field  response  to  a  point  dipole  current  source,  e  (r,  co)  is  the  dielectric  function  and  I  is 
the  unity  tensor. 

We  employ  FDTD  calculations  [48]  to  solve  numerically  for  the  Green’s  tensors  of  our 
structures.  The  integral  of  equation  (A.l)  is  evaluated  by  adapting  a  procedure  established 
in  [36]  and  by  using  a  deformed  contour  co (£)  =  £  J\  +icr/£  in  the  upper  half  of  the 
complex  frequency  plane,  parameterized  by  a  real  number  i;  J  0  and  a  constant  o  >  0.  As 
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explained  in  [36],  this  is  equivalent  to  solving  the  Green’s  tensor  at  real  frequencies  £  with 
a  fictitious  global  conductivity  applied  to  the  dielectric  function  e'(r,  £)  =  (1  +  kr/£)e(r,  £). 
The  integration  can  then  be  performed  in  the  time  domain  (via  the  convolution  theorem)  and 
converges  quickly  due  to  fast  decay  from  o . 

Specifically,  equation  (A.l)  is  numerically  evaluated  using 

h  C°° 

Uc p(ra)  =  —  /  dt  Im  [g/iV(— 0]  */r  •  ESC;V(ra,  t),  (A.l) 

2.K  Jo 

where  Esc  y(ra,  t)  is  the  (real)  electric  field  generated  by  a  point  dipole  current  source  J  =  8  (t)xv 
(xv=x,y,z)  located  at  the  position  ra  and  scattered  by  a  structure  with  a  dielectric  function 
e'(r,  £).12  Here,  the  indices  ji  and  v  are  repeated  for  summation  convention,  gllv(t)  is  the 

Fourier  transform  of  g,lv(£)  —  — i£^l  +  j(l  +  ]|)©(£)ar°v((z>(§))  and  ©(£)  is  the  Heaviside 
step  function. 

A.l.  Material  dispersion 

For  the  initial  calculations,  we  take  the  dielectric  constant  e  to  be  frequency  independent, 
e(r,  oj)  — »■  e(r).  In  principle  it  is  important  to  take  into  account  the  frequency  dependence  of 
the  refractive  index  when  calculating  the  CP  potential.  However,  for  the  trapping  of  atoms  at 
distances  of  lOOnm  from  dielectric  structures,  the  dominant  contribution  to  the  force  comes 
from  the  material  response  near  the  atomic  resonance  frequencies.  This  is  dictated  by  the  smooth 
frequency  cutoff  coming  from  the  atomic  polarizability;  see  equation  (A.l),  as  well  as  [49]  for 
more  detailed  discussions.  At  distances  on  the  order  of  the  resonant  wavelengths  or  more  (not 
relevant  to  our  calculations),  the  low-frequency  material  response  would  become  increasingly 
important  [49].  Since  the  refractive  index  of  silicon  nitride  is  approximately  constant  across 
the  optical  domain,  we  adopt  the  approximation  based  on  a  frequency  independent  dielectric 
constant  e(r). 

In  order  to  fully  justify  our  assumption,  we  have  performed  calculations  using  full  material 
dispersion,  and  have  found  the  results  well  approximated  by  calculations  using  constant  n  —  2. 
For  the  simulated  material  in  FDTD  calculations,  we  have  taken  the  measured  dispersion  in 
dc,  THz  [50]  and  optical  [51]  regimes  and  use  a  Lorentz  oscillator  model  to  fit  the  complex 
refractive  index  n  =n+ik\  see  figures  (A.la)-(b).  The  resulting  CP  potentials  are  shown  in 
figure  A.  1(c)  for  the  single  nanobeam  structure  and  in  figure  A.  1(d)  for  the  double  nanobeam 
structure.  Each  of  them  differs  from  the  values  calculated  with  n  =  2  by  no  more  than  8%  both 
at  the  trap  center  and  closer  to  the  walls. 

The  time  required  for  computation  increases  linearly  with  the  number  of  Lorentz  oscillators 
used  to  model  the  material,  and,  for  the  nine  oscillator  model  used  here,  the  computation  time 
increases  by  an  order  of  magnitude.  Due  to  limited  computational  resources,  we  have  used  the 
frequency  independent  results  as  a  reasonable  approximation  when  calculating  the  CP  potential 
throughout  both  structures. 

12  To  obtain  the  CP  potential  for  a  single  atom  inside  a  periodic  structure,  we  impose  Bloch-periodic  boundary 
conditions  on  a  single  unit-cell  in  the  FDTD  calculations.  We  then  sum  over  all  fields  ESCiV(ra,  t;  kx)  with  the 
Bloch  wavevector  kx  across  the  first  Brillouin  zone  to  obtain  ESCiV(ra,  t)  =  gg  /-i/ad£*ESCjV(ra,  t;kx),  which  is 
mathematically  equivalent  to  solving  for  the  fields  of  a  single  point  dipole  in  an  infinite  periodic  structure.  The 
integral  is  approximated  using  a  10  point  Gaussian  quadrature. 
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Figure  A.1.  Panels  (a)  and  (b)  show  the  real  («)  and  imaginary  ( k )  refractive 
indices,  respectively,  of  silicon  nitride  obtained  from  a  model  of  9  Lorentz 
oscillators.  Panel  (c)  shows  the  CP  potential  calculated  for  the  single  nanobeam 
structure  along  the  a -axis.  The  blue  line  shows  the  calculation  using  a  frequency 
independent  n  —  2.  At  the  discrete  points  x  —  —50,  0  and  50  nm,  the  frequency 
independent  results  (red  circles)  were  compared  with  results  using  the  full 
frequency  dependent  material  response  (black  diamonds).  Panel  (d)  shows  a 
similar  plot  along  the  y-axis  of  the  double  nanobeam  structure. 


Appendix  B.  Calculation  of  yto t 


To  determine  the  total  spontaneous  decay  rate  ytot  for  an  atom  in  our  structures,  we  also  solve 
for  the  classical  Green’s  tensors  and  evaluate  ytot  via  [34,  35,  46] 


2/To<r>y 

Ytot  =  —i — Im 
n 


5^Tr[Dr  G(ra,  ra,  coj)] 
{0) 


(B.l) 


where  D;  =  ( { 0}  | d:|  y )  ( y  | d |  {0} )  is  the  dipole  matrix  element  between  the  ground  state  manifold 
and  the  excited  state  j,  and  coj  is  the  transition  frequency.  The  total  decay  rate  ytot  —  yiD  +  y' 
includes  the  decay  rate  ylD  to  a  guided  mode  of  interest  as  well  as  the  rate  y'  to  all  other  modes 
of  the  structure,  including  lossy  modes.  As  discussed  below,  the  contributions  of  yjD,  y'  to  ytot 
can  be  estimated  from  the  global  frequency  dependence  ytot(oj). 

To  obtain  ytot((o),  we  evaluate  the  Green’s  tensor  for  the  real  dielectric  function  e(r)  using 
the  FDTD  method,  followed  by  a  discrete  Fourier  analysis. 


Appendix  C.  Validation 

To  validate  our  numerical  procedures,  we  have  performed  calculations  of  Uc p  for  several 
geometries  where  analytical  solutions  are  available,  including  an  atom  near  an  infinite  dielectric 
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or  metal  half-space  [52]  and  an  atom  located  above  an  infinite  dielectric  grating  [53],  and  found 
excellent  agreement  between  our  simulations  and  the  exact  results. 

We  have  validated  our  calculations  of  ytot  for  the  cases  of  an  (atomic)  dipole  near 
an  infinite  dielectric,  metallic  parallel  plates,  a  nanofiber  [54]  and  2D-photonic  band-gap 
microcavities  [55]. 

Appendix  D.  Guided  mode  resonances 

For  our  structure  with  an  infinite  number  of  unit  cells,  a  guided  mode  (denoted  by  A)  contribution 
to  the  imaginary  part  of  the  Green’s  tensor  can  be  calculated  as  [46,  56],  Im  [GjD(ra,  ra,  <y)]  = 
£/c'2U;  (ra;  kx)  ®  u*(ra;  kx)/2covx,  when  the  frequency  co  intersects  the  frequency  band  cok  at  a 
propagation  constant  kx  below  the  light  line.  Here,  u>  (r;  kx)  is  the  orthonormal  mode  function, 
and  vk  is  the  group  velocity,  both  available  via  numerical  calculations  [44]  (see  also  footnote  8). 
As  we  scan  the  frequency  co,  Im  [GjD]  increases  monotonically  and  diverges  as  kx  approaches 
the  A -point,  where  vk  — »•  0.  The  guided  mode  Green’s  tensor  vanishes  when  co  lies  beyond  the 
frequency  of  the  band  edge. 

Based  on  this  analysis,  we  can  evaluate  the  decay  rate  yip*  (co)  into  the  designated  probe 
mode  for  an  infinite  structure,  and  compare  it  with  the  heights  of  resonant  features  in  yto t(co)  for 
finite  structures  with  different  numbers  of  unit  cells,  as  shown  in  figures  4(b)  and  5(b).  Indeed, 
the  actual  yiD  of  a  finite-size  structure  must  deviate  from  y[^}  due  to  boundary  conditions 
that  transform  a  continuous  spectrum  into  a  discrete  set  of  resonant  peaks  [40],  as  shown  in 
figures  4(b)  and  5(b).  When  the  number  of  unit  cells  is  increased  in  our  calculation  of  ylD  over 
the  range  N  —  1 1  to  81,  we  find  that  the  frequencies  co(n)  of  the  resonant  peaks  shift  in  position 
and  the  peaks  change  height.  As  documented  by  the  black  diamonds  in  figures  1(b),  (f),  the  co(n) 
arise  from  the  discrete  set  of  propagation  constants  kf 1  ~  itn/aN  imposed  by  the  boundary 
conditions  for  the  finite  structures  with  n  either  even  or  odd. 

The  peaks  in  ytot  at  the  set  of  frequencies  coin}  build  up  on  top  of  a  fairly  constant 
background  within  the  frequency  range  displayed  in  figures  4(b)  and  5(b).  We  assume  that  this 
background  represents  the  contribution  of  y'  to  ytot,  and  subtract  the  background  to  estimate 
yiD-  The  resulting  form  for  y\o(co)  consists  of  a  set  of  resonant  peaks  whose  heights  at  discrete 
coM  qualitatively  map  out  y1(^o)  calculated  for  the  infinite  structure,  with  the  maximum  peak 
height  for  ylD  occurring  for  the  peak  closest  to  the  band  edge.  Moving  further  away  from  the 
band  edge,  we  find  that  our  numerical  estimate  of  ytot(co)  —  y'(co )  asymptotes  to  the  calculated 
value  of  y1(^o)  (co)  reasonably  well. 

D.l.  Estimation  of  yt  D  and  y' 

From  the  previous  discussions,  we  identify  that  the  decay  rate  into  other  modes  y'  can  be  read 
off  from  the  broad  background  in  ytot.  Specifically,  we  estimate  y'  —  yVA(oj)  at  a  frequency  co’ 
just  across  the  band  edge  and  away  from  any  resonant  peak  for  a  guided  mode.  The  decay  rate 
into  the  probe  mode  yiD  can  then  be  estimated  using  y1D  =  ytot  —  y'. 

For  the  single-beam  structure  and  the  atomic  spin  orientation  shown  in  figure  4(b),  we 
find  a  peak  total  decay  rate  ytot/yo  ~  15  and  a  background  level  y'/yo  ^  1.2  near  the  cesium 
D2-line  frequency  vA  =  352  THz.  We  estimate  the  coupling  to  the  resonant  probe  mode 
ylD  —  ytot  —  y '  14yo-  For  the  double-beam  structure  and  spin  orientation  shown  in  figure  5(b), 
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we  find  Ytot/Yo  ~  22,  y'/yo  ~  1,  and,  therefore,  ym/Yo  ~  21.  Here,  y^/ln  =  5.2MHz  is  the 
free-space  (vacuum)  decay  rate  for  the  D2  line. 

Appendix  E.  Effective  area  and  mode  volume  for  probe 

Both  the  single  and  double  nanobeam  structures  lead  to  atom  localization  in  a  region  of  large 
amplitude  for  the  probe  field.  One  measure  of  the  strength  of  the  atom-field  coupling  is  the 
effective  mode  volume  Vm  per  unit  cell,  where 

Vm  =  J  €(r)|E(r)|2d3r/e(rmin)|E(rmin)|2.  (E.l) 

Here  the  integration  is  carried  out  over  the  volume  of  a  unit  cell.  That  is,  the  integration  domain 
along  propagation  direction  x  extends  over  the  distance  a  (i.e.  the  lattice  constant),  while  in  the 
transverse  y,  z  directions,  the  integration  domain  is  from  —  oo  to  +00. 

For  the  single  nanobeam,  the  probe  mode  has  a  global  maximum  at  rmin  =  0  and  an  effective 
mode  volume  Vm  ~  0.13  gm3.  For  the  double  nanobeams,  the  probe  mode  has  a  saddle-like 
intensity  distribution  around  rmjn  =  0,  resulting  in  Vm  ~  0. 1 1  gm3  for  a  unit  cell. 
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